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ABSTRACT 

A new tool for the Evolutionary Synthesis of Stellar Populations (EPS) is presented, 
which is based on three independent matrices, giving respectively: 1) the fuel con- 
sumption during each evolutionary phase as a function of stellar mass; 2) the typical 
temperatures and gravities during such phases; 3) colours and bolometric corrections 
as functions of gravity and temperature. The modular structure of the code allows to 
easily assess the impact on the synthetic spectral energy distribution of the various 
assumptions and model ingredients, such as, for example, uncertainties in stellar evolu- 
tionary models, mixing length, the temperature distribution of horizontal branch (HB) 
stars, AGB mass loss, and colour-temperature transformations. The so-called "AGB- 
Phase Transition" in Magellanic Cloud clusters is used to calibrate the contribution of 
the Thermally Pulsing Asymptotic Giant Branch phase (TP- AGB) to the synthetic in- 
tegrated luminosity. As an illustrative example, solar metallicity (Y — 0.27, Z = 0.02) 
models, with ages ranging between 30 Myr and 15 Gyr and various choices for the 
slope of the Initial Mass Function (IMF), are presented. Synthetic broad band colours 
and the luminosity contributions of the various evolutionary stages are compared 
with LMC and Galactic globular cluster data. In all these cases, a good agreement is 
found. Finally, we show the evolution of stellar mass-to-light ratios in the bolomet- 
ric and U,B,V,R, and K passbands, in which the contribution of stellar remnants is 
accounted for. 



Key words: 

AGB 



galaxies: star clusters: evolution - galaxies: Magellanic Clouds - stars: 



1 INTRODUCTION 

In order to study the formation and the subsequent evolu- 
tion of galaxies, we have to understand the overall proper- 
ties of their stellar content. In a galaxy, a mixture of stellar 
populations of different ages and chemical compositions are 
present and many efforts have been made in recent years, 
inspired by the pioneer work of Tinsley (1980), aimed at 
modelling the spectral evolution of the various morphologi- 
cal types, from ellipticals to spirals (e.g. Bruzual 1983; Ari- 
moto & Yoshii 1987; Guiderdoni & Rocca-Volmerange 1987; 
Bruzual & Chariot 1993; Bressan, Chiosi & Fagotto 1994; 
Tantalo et al. 1996). 

The building blocks of Evolutionary Population Synthe- 
sis (EPS) are models for Simple Stellar Populations (SSPs), 
that are assemblies of chemically homogeneous and coeval 
single stars. Thus, before facing the problem of modelling 
a complex stellar population, as a galaxy, the preliminary 
condition is to check the accuracy of SSPs computations 
and their adequacy to reproduce the observable features of 
stellar clusters, that better resembling the definition of SSP. 



Two main approaches have been followed in computing EPS 
for SSPs. The 'Isochrone Synthesis' technique (Chariot & 
Bruzual 1991) consists in summing up the contributions to 
the flux in the various passbands of all mass-bins along one 
isochrone, after assuming an Initial Mass Function (IMF). 
The integration starts from a lower mass limit and ends at 
the latest mass point on the isochrone itself, that usually co- 
incides with the end of the so-called Early Asymptotic Giant 
Branch (E-AGB) phase. Later stellar phases are then added 
following individual recipes (e.g. Chariot & Bruzual 1991). 
An alternative approach to compute the luminosity contri- 
butions of Post-Main Sequence evolutionary stages is based 
on the so-called Fuel Consumption theorem (FCT, Renzini 
& Buzzoni 1986). Here the main ingredient of synthesis is 
the amount of nuclear fuel (i.e. the hydrogen and/or helium 
mass) that is burned in each evolutionary stage. The fuel, 
a natural product of evolutionary stellar sequences, is con- 
verted into bolometric luminosity thanks to the FCT. The 
luminosities in the various passbands are then computed 
thanks to a temperature/gravity-colour set of transforma- 
tions. In principle the two approaches are equivalent, as both 
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are based on stellar evolutionary models, but their practi- 
cal implementations may result in sizable differences. In the 
isochrone synthesis method, the integration variable is the 
stellar mass, that varies along the isochrone. So, bright, but 
short-lived evolutionary stages span a very narrow range 
in stellar mass and integrations must use very small AM 
steps. (We have noticed that mass steps as small as 1CF 5 
Mq are required to achieve a few per cent stability in nu- 
merical results.) Coarse mass zoning and poor interpolation 
algorithms are responsible for the sizable fluctuations oc- 
casionally exhibited by population models constructed with 
the isochrone synthesis method (see e.g. Fig. 7 in Chariot 
& Bruzual 1991). The method based on the FCT is instead 
much more robust in this respect, as the integration vari- 
able, the fuel, is directly proportional to the contribution of 
the various phases to the total luminosity. 

At present this method has been followed by Buzzoni 
(1989) to construct SSPs for ages older than 4 Gyr. In this 
work we extend the approach to young and intermediate age 
populations and present SSP models for solar metallicity and 
ages from 30 Myr to 15 Gyr. The synthesis computational 
code has been constructed with the aim of better controlling 
the influence of various model ingredients on SSP models. 
The comparison with Magellanic Clouds and Galactic glob- 
ular clusters allows a detailed test for the results. The paper 
is organised as follows. In Section 2 a detailed description 
of the algorithm is presented; Section 3 illustrates all model 
ingredients: stellar input and temperature-colour transfor- 
mations. Section 4 describes the results and Section 5 shows 
the various observational tests performed. Finally, in Section 
6, our conclusions are drawn. 



2 THE ALGORITHM 

The EPS code used in this paper is based on three inde- 
pendent sets of matrices containing the ingredients of the 
synthesis, namely: 

(i) The energetics: the nuclear fuel burned during each 
evolutionary phase as a function of stellar mass; 

(ii) The surface parameters: i.e. the effective tempera- 
tures and surface gravities during such phases; 

(iii) The transformations to observables: i.e. colours and 
bolometric corrections as functions of gravity and tempera- 
ture. 

Following Renzini & Buzzoni (1986; hereafter RB86), 
the total bolometric luminosity of a Simple Stellar Popula- 
tion (SSP) of age t can be written as the sum of two terms, 
one for the core hydrogen burning stars (the Main Sequence 
(MS) stars), one for the evolved ones (the Post-Main Se- 
quence (PMS) stars) 



L T ° (t) — Lms(£) + £pMsM- 



(1) 



Lms depends on the adopted mass-luminosity relation 
L(M,t) and Initial Mass Function (IMF) *(M) = 
AM~^ 1+X ^ (A being the scale factor giving the size of the 
stellar population). Lpms is directly connected with the total 
nuclear fuel burned by stars populating the various evolved 
stages, through the FCT (see 2.1). 



To evaluate the MS bolometric luminosity, we need a set 
of isochrones, on which to perform a numerical integration 
weigthed on the IMF 



r bol 



/■M TO (t) 

(t) = / L(M,i)#(M)dM 

J M int 



(2) 



The lower integration limit is the minimum mass igniting H 
in the core (Mi n f). Mto (t) (the turnoff mass) is the mass in 
the verge of exhausting hydrogen in the core and leaving the 
MS. We chose as Mto (t) the model with the maximum ef- 
fective temperature along the isochrone as, by definition, the 
turnoff point is the bluest point on the Hertzprung-Russel 
(H-R) diagram. The relation between the age and the corre- 
sponding turnoff mass is the evolutionary clock of a stellar 
population. 



2.1 Post-Main Sequence contributions: the Fuel 
Consumption Theorem 

The Fuel Consumption theorem (FCT) states that: the con- 
tribution of stars in any given Post-MS stage to the inte- 
grated bolometric luminosity of a SSP is directly proportional 
to the amount of fuel burned during that stage (RB86). Its 
analytical form is 



L)°\t) = 9.75 x lQ w b{t)F s {M TO ) 



(3) 



where L^ ol (t) is the total bolometric luminosity of stars pop- 
ulating the PMS stage j at the age t. b(t), the "evolutionary 
flua?' , is the number of stars evolving off the Main Sequence 
per year. fj(Mro) is the nuclear fuel (in Mq) burned by 
stars with M=Mto in their j-th evolutionary PMS phase. 
The function b(t) is defined as 



b(t) = *(M TO )| 



rfMxo 
dt 



(stars /yr) 



(4) 



where *(M T o) is the IMF computed for M = Mto and 
Mto is the time derivative of the relation Mxo(i)- At any 
SSP age, the mass of evolved stars is only slightly greater 
than Mto(£)- Hence, to a fair approximation b(t) also gives 
the rate at which stars enter or leave any specific post-MS 
evolutionary stage. 

The fuel consumption F- 3 is given by 

Fj =m j H + 0.1m j Ho (Mq) (5) 

where mj H and mj Ho are the amount of hydrogen and he- 
lium burned during the phase j, the coefficient 0.1 taking 
into account that the energy release from helium-burning is 
~ 1/10 than from H-burning. The close coincidence in mass 
between evolved stars and Mro(t) allows us to refer the 
PMS fuels to the values proper to a mass equal to Afro(t)- 
As a consequence, mass and age indices are equivalent for 
PMS stars and the index j will be used to label star masses 
in the MS and evolutionary stages in the PMS. 

This approximation, in which one assumes the flux b(t) 
to be constant through the various PMS stages, was exten- 
sively criticized by Chariot & Bruzual 1991. However, it has 
been demonstrated that this is actually a very accurate ap- 
proximation (Renzini 1994). 
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2.2 From Bolometric to Monochromatic 



To compute the SSP total luminosity in the generic photo- 
metric band A, i.e. 



L^(t)=L x M3 {t)+L x PM3 {t), 



(6) 



we need temperature/gravity - bolometric correction/colour 
sets of transformations. The (T e ,g)ij matrix contains the 
effective temperatures and the surface gravities of stars in 
the j-th phase, at the SSP age i (see 3.3 for a detailed de- 
scription of adopted values). The colour index/bolometric 
correction matrix (CI, BC)i m contains the trasformations 
to the observables, i.e. bolometric corrections and colours 
as functions of the surface parameters (T e ,g). I is an effec- 
tive temperature index and m refers to surface gravity: we 
defer to 3.4 for a detailed description of adopted sets. The 
A-luminosity of the j'-star at the SSP age i, L^, is a fraction 
of its bolometric luminosity L^" 1 , as 

- /-!f 04) (7) 

Luminosities in eq. (7) are in solar units and Qy is defined 

as 



„\ -0.4f.BCi.- 

Q x = 10 1 



(8) 



where BC-^ is the bolometric correction to the A-band and 
BCq the same for the Sun. 

In order to get the luminosity contributions of PMS 
phases, we divide each phase j into a certain number of 
photometric sub-phases, burning an equal fraction of the 
total fuel Fij and having a constant Q-factor (eq. 8). Then 
we apply eq. (7), summing over these sub-phases. The SSP 
total A-luminosity at the age i is finally 



(9) 



A correct link between MS and PMS luminosity contribu- 
tions is ensured by adopting the same normalization con- 
stant A for *(M) in eq. (2) and (3). 

The modular structure of the code allows to easily 
change a particular ingredient and/or assumption, indepen- 
dently of the others, as assessing just its own impact in the 
final output (see 5.4 for some examples). 



3 MODEL INGREDIENTS 

In this section we describe all adopted model ingredients. 
The stellar stages included in the synthesis range from the 
Main Sequence to the end of the Thermally Pulsing Asymp- 
totic Giant Branch (TP-AGB). The Post-AGB has not been 
taken into account as the UV output of old populations is 
extremely model dependent (cf. Greggio & Renzini 1990). 



3.1 Stellar input 

For the MS mass-luminosity relation L(M,t) (eq. 2), a ho- 
mogeneous set of solar chemical composition (Y — 0.27, Z = 
0.02) isochrones (Castellani, Chieffi & Straniero 1992, here- 
after CCS92), with ages ranging from 30 Myr to 10 Gyr has 
been adopted. An additional isochrone for 15 Gyr has been 
constructed from the same set of stellar tracks kindly pro- 
vided by Oscar Straniero. The dwarf-MS component down 
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Figure 1. The relation Mxo(*) an d its time derivative Mxo (t in 
yr, Mxo i n Mq). The relation has been obtained via an analytical 
fit on CCS92 isochrones. 



to 0.1 Mq, has been added using the models of VandenBerg 
et al. (1983) for Y=0.25, Z=0.02. The small difference in 
the helium content is not crucial in this context, as the loca- 
tion of the lower Main-Sequence, in the logi - logT c plane, 
strongly depends on metallicity, but it is rather insensitive 
to the small variations of Y parameter. 

The relation Mxo /age (Fig. |l|) is obtained 
lytical fit on CCS92 isochrones 



via an ana- 



log -Mxo (t) = alogt + blog 2 t + clog 3 t + d 



(10) 



where a = 0.212567, b = - 0.108394, c = 0.005737, d = 
2.981730 (M TO in Mq, t in yr). The time derivative of 
Mxo(i) is also displayed in Fig. hi. As it is evident, the SSP 
evolutionary clock, Mxo, is very fast for the large values of 
Mxo, then its rate drops greatly as Mxo decreases. In fact, 
M TO ranges from ~ 8.3 at t = 30 Myr to ~ 0.94 at * = 15 
Gyr while, in the same age range, Mxo drops by more than 
4 orders of magnitude. This implies a faster SSP luminosity 
evolution at young ages (t ^ 1 Gyr). 

The stellar evolutionary tracks used to derive eq. (10) 
have been constructed using canonical assumptions con- 
cerning convective boundaries, i.e. no overshooting. Slightly 
higher values of Mto for given age are obtained when allow- 
ing for convective overshooting, the size of the effect depend- 
ing on a free parameters. However, no general agreement has 
yet emerged as to what value should be assigned to the free 
parameter, and therefore we prefer to use canonical models. 

Fig. |2| shows the evolutionary flux b(t), i.e. the rate of 
evolution off the Main Sequence (eq. 4), for four different 
choices of the IMF slope. The functions b(t) are normalized 
at b — 1 for t — 15 Gyr, thus a comparison between the dif- 
ferent slopes is possible regardless of the SSP size. b(t) is the 
product of two terms: one increasing (ty) and one decreasing 
(Mxo) with age. The flatter the IMF, the slower the increase 
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Figure 2. The evolutionary flux b(t) for three choices of the IMF 
slope x. The effect of a multislope, like a Scalo type, is also shown, 
in which x is -1 for M < 0.3, 1.35 for 0.3 < M < 2, 1.7 otherwise. 
Functions b(t) are normalized at t = 15 Gyr. 

in time of <&(M) and the rapid decrease in Mto dominates. 
This results in a greater SSP luminosity evolution for flatter 
IMFs. 

3.2 The Fuel Consumption matrix Fy 

Table 1 is our Fy matrix. It gives the amount of hydrogen 
and helium (in Mq) that are burned by a star of mass Mxo,i 
in each of its PMS stages. The total nuclear fuel follows from 
eq. (5), while JWro.i and the corresponding age are related 
by eq. (10). The adopted values of Fy in Table 1 are basi- 
cally the same as those used by Renzini (1992). They have 
been derived from various sets of evolutionary sequences, in- 
cluding Sweigart & Gross (1978), Becker (1981), Sweigart, 
Greggio & Renzini (1989), Renzini & Voli (1981), and Greg- 
gio & Renzini (1990). These data have been complemented 
with data from Chieffi & Straniero (1992, private commu- 
nication) for the fuel consumption during the RGB phase. 
In principle it would have been most appropriate to use fuel 
consumption data from a unique, homogeneous set of evo- 
lutionary calculations. Unfortunately, such a complete set, 
covering all evolutionary phases for all stellar masses is still 
not available. However, the present compromise is largely 
adequate for the scope of this paper. 

A special care was devoted to estimate the fuel con- 
sumption during the TP-AGB phase. Uncertainties in mass 
loss, mixing, and efficiency of hydrogen burning at the bot- 
tom of the convective envelope (known as envelope burning 
process) together prevent pure theory from predicting the 
amount of fuel burned during this stage (RB86). The uncal- 
ibrated models of Renzini and Voli (1981, hereafter RV81) 
used by RB86 largely overestimate the TP-AGB contribu- 
tion for populations younger than ~ 10 s yr in comparison 
with the observed contribution in Magellanic Cloud clusters 
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Figure 3. The synthetic AGB bolometric contribution compared 
with LMC GCs data from FMB90. E-AGB and TP-AGB con- 
tributions are also shown separately. Each data point has been 
obtained through a procedure aimed at minimizing the stochasti- 
cal fluctuations between clusters having an equal SWB parameter 
(see the text). Error bars indicate the r.m.s. 



in the corresponding age range (Frogel, Mould & Blanco 
1990, hereafter FMB90). 

The origin of the discrepancy has been controversial for 
several years, but a decisive hint came from the discovery 
that TP-AGB stars experiencing the envelope burning pro- 
cess do not follow the canonical AGB core mass-luminosity 
relation (Blocker & Schonberner 1991). The 7 Mq model by 
Blocker & Schonberner is much more luminous and evolves 
~ 10 times faster than that predicted by the standard core 
mass-luminosity relation (e.g. Paczynski 1970a, b) that was 
used by RV81. The consequence is a drastic decrease of 
the TP-AGB lifetime and of the TP-AGB fuel consump- 
tion, that results in a lower contribution of this phase to the 
integrated light of a population. 

The transition between TP-AGB stars that experience 
the envelope burning process and those which do not is very 
sharp, and takes place at a threshold mass that depends 
sensitively to the mixing length parameter a (RV81). Now, 
Magellanic Cloud clusters offer an opportunity to calibrate 
this parameter in such a way to reproduce the observed con- 
tribution of AGB stars. Following Blocker & Schonberner 
(1991) we assume that TP-AGB stars experiencing enve- 
lope burning consume ~ 10 times less fuel than predicted 
by RV81, while we retain RV81 fuel consumptions for stars 
that do not experience envelope burning. In this way the cal- 
ibration procedure reduces to find the value of a that best 
accounts for the data. The result is shown in Fig. ^ which 
compares the theoretical AGB contributions for a — 2 to 
the Magellanic Cloud cluster data from FMB90. For this 
value of a stars more massive than ~ 3 Mq experience the 
envelope burning process and have the TP-AGB fuel con- 
sumption correspondingly reduced. 
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Table 1. The Fij matrix: the masses of hydrogen (mjj) and helium (mjje) (both in Mq) burned in 
each evolutionary PMS j phase, namely SGB, RGB, HB, E-AGB, TP-AGB, by a star of M = M TO ,i 
(i is the age index). Values are the same as in Renzini (1992; see the text). 



SGB RGB HB E-AGB TP-AGB 

AGE (Gyr) Mxo i fnB ma "m-H f"He ma friHe m H mn c 



0.03 


8.2247 


0.0222 


0.0058 


0.3116 


0.8280 


0.0000 


0.2091 


0.0339 


0.0491 


0.06 


5.9188 


0.0267 


0.0059 


0.1898 


0.5211 


0.0000 


0.4210 


0.0339 


0.0491 


0.1 


4.7334 


0.0296 


0.0046 


0.1482 


0.3993 


0.0015 


0.4567 


0.0339 


0.0491 


0.2 


3.6001 


0.0323 


0.0031 


0.1233 


0.3115 


0.0081 


0.3918 


0.0339 


0.0491 


0.4 


2.7777 


0.0331 


0.0041 


0.1109 


0.2647 


0.0110 


0.3448 


0.2071 


0.3138 


0.6 


2.3977 


0.0329 


0.0108 


0.1136 


0.2366 


0.0181 


0.2971 


0.1931 


0.2927 


0.8 


2.1640 


0.0315 


0.0388 


0.1108 


0.2201 


0.0153 


0.2867 


0.1744 


0.2643 


1 


1.9985 


0.0307 


0.0993 


0.0815 


0.2349 


0.0119 


0.2802 


0.1427 


0.2163 


1.5 


1.7408 


0.0364 


0.1453 


0.0668 


0.2386 


0.0107 


0.2969 


0.0979 


0.1483 


2 


1.5784 


0.0431 


0.1595 


0.0621 


0.2373 


0.0113 


0.3083 


0.0747 


0.1132 


3 


1.3844 


0.0509 


0.1689 


0.0646 


0.2472 


0.0139 


0.2921 


0.0463 


0.0701 


4 


1.2651 


0.0572 


0.1806 


0.0618 


0.2594 


0.0116 


0.2827 


0.0363 


0.0550 


6 


1.1341 


0.0625 


0.1885 


0.0595 


0.2665 


0.0090 


0.2750 


0.0254 


0.0384 


8 


1.0588 


0.0653 


0.1921 


0.0583 


0.2693 


0.0076 


0.2710 


0.0191 


0.0289 


10 


1.0110 


0.0660 


0.1947 


0.0560 


0.2671 


0.0086 


0.2671 


0.0151 


0.0228 


15 


0.9380 


0.0659 


0.1991 


0.0508 


0.2598 


0.0124 


0.2598 


0.0090 


0.0136 



In Magellanic Cloud clusters the number and luminos- 
ity of AGB stars is subject to large fluctuations from one 
cluster to another of similar age, and therefore so does the 
AGB contribution to the total light. This is a result of small 
number statistics, as clusters contain at most just a few AGB 
stars. To cope with this problem, clusters have been grouped 
in age bins following the classification of Searle, Wilkinson 
and Bagnuolo (1980), with the SWB type being related to 
age by the relation given by FMB90. Within each bin the 
AGB contribution to the bolometric luminosity was then 
obtained by adding together the luminosity of all the AGB 
stars in the binned clusters, and dividing by the sum of all 
the integrated luminosities of the same clusters. Fig. [] shows 
that a jump in the AGB contribution takes place between 
SWB type III and SWB type V, or at an age of ~ 0.2 Gyr, 
with the AGB contribution suddenly increasing from ~ 10 
per cent to over ~ 30 per cent. This corresponds to the 
so-called "AGB phase transition" (RB86) due to the devel- 
opment of a bright and well populated TP-AGB phase (see 
also Renzini 1992). 

The TP-AGB phase is populated by C and M spectral 
type stars. The production of C (carbon) stars is a func- 
tion of both metallicity and the initial stellar mass, hence of 
the SSP age (cf. RV81). From FMB90 data for LMC clus- 
ters (their Table 4) we have evaluated the C and M star 
luminosity contributions to the total light for each age bin. 
As 

^tp-agb -Ftp -agb (eq. 3), we have then calibrated 
the fractions of TP-AGB fuel burned by M and C stars as 
functions of SSP age. The use of LMC GCs as calibrators 
is unavoidable, as they represents the only available sample 
of clusters covering a broad age-range and being populous 
enough, so ensuring that the AGB is reasonably well sam- 
pled. These clusters have typically ~ 1/2 Zq metallicity, 
but a similar sample for solar metallicity is not available. 
The dependence on Z may be taken into account using the 
prescriptions of RV81. This will be discussed in sec. 5.4.2. 

The calibration is shown in Table 2 (and Fig. Q). As it 
can be noticed, carbon stars characterize intermediate age 



Table 2. The calibration of the TP-AGB fuel (see the text). 
Column 2 give the total TP-AGB fuel, Column 3 and 4 its fraction 
burned by C and M type stars respectively, as functions of the 
SSP age. 



AGE(Gyr) 


-Ftp -agb 


%-Ftp-agb 


%F M 

'"^TP-AGB 


0.03 


0.0388 


0.0000 


1.0000 


0.06 


0.0388 


0.0000 


1.0000 


0.1 


0.0388 


0.0000 


1.0000 


0.2 


0.0388 


0.0000 


1.0000 


0.4 


0.2385 


0.1138 


0.8862 


0.6 


0.2224 


0.1138 


0.8917 


0.8 


0.2008 


0.1160 


0.9019 


1 


0.1643 


0.1421 


0.8579 


1.5 


0.1127 


0.1065 


0.8935 


2 


0.0860 


0.0846 


0.9154 


3 


0.0533 


0.0689 


0.9311 


1 


0.0418 


0.0583 


0.9417 


6 


0.0292 


0.0657 


0.9343 


8 


0.0220 


0.0000 


1.0000 


10 


0.0174 


0.0000 


1.0000 


15 


0.0104 


0.0000 


1.0000 



SSPs, those having a developed TP-AGB. Very old stellar 
populations (e.g. galactic globular clusters and the Galactic 
Bulge) do not have C stars, in agreement with theoretical 
expectations. 

Fig. |H| shows the total amount of fuel burned in the 
whole PMS, Fj.T = -Fij, as a function of the SSP age and 
the turnoff mass. The total fuel decreases with the decreas- 
ing mass, till the onset of AGB phase transition, at ages 
around 0.2 Gyr. Starting from t <; 1 Gyr, Fi,T tends to re- 
main nearly constant as the decrease in the TP-AGB fuel is 
partly compensed by the increase in the fuel burned on the 
red giant branch (RGB). This epoch is referred to as "RGB 
phase transition" (RB86), the mass at the transition being 
around ~ 2.5 Mq, in dependence of the chemical composi- 
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Figure 4. The total amount of TP-AGB fuel (solid line), its 
fraction burned by C type stars (dotted line ) and by M type 
stars (dashed line ), as functions of age (see the text). 



Figure 6. The fuel consumption in each evolutionary PMS phase 
j as function of the SSP age; in panel (d), i*E— AGB (dotted line) 
and i^TP-AGB (dashed line) are separately shown. 
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increase of the RGB fuel consumption at t ~ 1 Gyr, which 
is due to the appearence of stars developing a degenerate 
helium core, a phenomenon known as RGB phase transition 
(RB86, Sweigart, Greggio & Renzini 1990). The small hump 
in the HB fuel consumption at nearly this same age is also 
due to the RGB phase transition, with stars with a minimum 
value of the core mass populating the HB, hence experienc- 
ing a prolongued HB phase leading to a local maximum in 
the hydrogen consumption (cf. Table 1). 

Finally, Fig. ^ shows the contribution of various evolu- 
tionary stages to the total bolometric luminosity, Ly /L^ t, 
obtained with the procedure described in sec. 2 and the in- 
gredients of sec. 3.1 and 3.2. The energetic of young SSPs 
(t ;SS 0.3 Gyr) is dominated by MS stars, that of intermedi- 
ate age ones (0.3 Ss t SS> 3 Gyr) by AGB stars and, finally, by 
RGB stars for t J; 3 Gyr. For ages older than ~ 6 Gyr, the 
relative MS contribution tends to rise. This effect is caused 
by the MS integrated luminosity decreasing slower than b(t). 
As the total fuel keeps nearly constant, the net result is a 
lower total luminosity as mainly due to a lower PMS lumi- 
nosity. 



Figure 5. The total fuel consumption (in Mq) during the whole 
post-MS, as a function of the SSP age (lower scale) and the cor- 
responding turnoff mass (upper scale). 



tion. As it can be noticed, the PMS total fuel depends only 
weakly on age, ranging from ~ 0.5 Mq at t ~ 30 Myr to 
~ 0.4 Mq at t ~ 15 Gyr. Also, it is possible to show that 
the MS luminosity contribution tends to be nearly constant 
after an age of ~ 1 Gyr. As a consequence, the luminosity 
evolution of a SSP is mainly controlled by the evolutionary 
flux b(t). Fig. | shows the time evolution of the fuel con- 
sumptions in the various phases. Worth noting is the sharp 



3.3 The temperature/gravity matrix (T e ,g)ij 

The (T c ,g)ij matrix contains the distribution of effec- 
tive temperatures and surface gravities of the evolutionary 
phases j at the SSP age i. 

Stellar effective temperatures depend crucially on the 
efficiency of convective energy transfer, parametrized by the 
mixing lenght parameter a. The use of uncalibrated the- 
oretical effective temperatures in EPS computations is ex- 
tremely dangerous, as the calibration of a is the synchroniza- 
tion of the EPS clock. CCS92 isochrones are computed for 
a = 1.6, a value that matches the observed Sun. Moreover, 
these isochrones have been tested on a sample of well-studied 
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Table 3. The (T e ,g);j matrix for the 30 Myr, the 1 Gyr and the 15 Gyr SSP: effective 
temperatures and surface gravities for the j evolutionary phases, as functions of the i ages. 
Values have been drawn from CCS92 isochrones. 
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Galactic clusters and account well for the observed loca- 
tion of stars in the various evolutionary stages (see CCS92) . 
The (T e , g)ij values for MS to E-AGB have been taken from 
CCS92 calibrated isochrones. V83 models have been used 
for the dwarf-Main Sequence component (0.1 < M < 0.6). 

For the TP-AGB phase, the relative proportions of C 
and M stars as functions of age (from FMB90, see sec. 3.2) 



are used as indicators of the evolution of the spectral type 
along the TP-AGB. 

The time evolution of the T c ranges spanned by PMS 
evolutionary stages is displayed in Fig. g, while Fig. |j shows 
the same for the MS. To note in Fig. [8| the red HB clump at 
late ages (see 5.4.1). As an example, Table 3 is the (T e ,gi)ij 
matrix for i=30 Myr, 1 Gyr and 15 Gyr (Mro = 8.2247, 
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Table 3 - continued 
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Table 4. The adopted solar normalization 

BCvq M bo i e M UQ M B q M V(S M R(S M KQ 
-0.08 4.75 5.61 5.48 4.83 4.31 3.41 

1.9985 and 0.9380 respectively). By changing the content of 
the (T c ,g) matrix, it is possible to infer the impact on the 
integrated photometric features of SSPs, of some model pa- 
rameters such as mixing length and mass loss, that directly 
influence effective temperatures and surface gravities (see 
5.4). 



(1966) empirical tables, on which a linear interpolation in 
(B — V) has been performed. For M type stars, the theoret- 
ical models by Bessel et al. (1989) have been used. Colours 
have been linearly interpolated in T c . For M and C type 
TP-AGB stars (V — R) and (V — K) values have been taken 
from various sets of observational data (Cohen et al. 1981; 
Cohen 1982; Westerlund et al. 1991). Finally, Table 4 shows 
the adopted values for the solar normalization (from Allen 
1991). By changing the {CI, BC)i m matrix, it is possible to 
explore the influence on model output of the photometric 
ingredients and the associated uncertainties. 



3.4 The colour /bolometric correction matrix 

The (CI, BC)\ m matrix contains bolometric corrections and 
colour indices, to distribute the total energy through the 
various passbands. It covers the range of (T e ,g)ij matrix 
and contains empirical colours suitable to TP-AGB M and 
C type stars. BCv and (17 — B) are taken from Kurucz 
(1992) theoretical grids, on which a quadratical interpola- 
tion in logTe, logy has been performed. For T a < 3500 °K, 
the limit of Kurucz models, some extrapolations have been 
necessary. (B — V) colours have been taken from CCS92 
isochrones. They are computed on the basis of the empirical 
color-temperature relation obtained by Arribas & Martinez 
(1988, 1989), for 4000 K < T e < 8000 K (see the authors 
for more details); outside this T c range, they are based on 
Kurucz models. (V — R) and (V — K) colours, for spec- 
tral types earlier than K, have been taken from Johnson 



4 RESULTS 

Evolutionary synthesis models for SSPs have been computed 
with the algorithm presented in Section 2 and the ingredi- 
ents described in Section 3. Models are characterized by a 
solar chemical composition and ages ranging from 30 Myr 
to 15 Gyr. The model predicts the SSP photometric proper- 
ties as functions of age, i.e. the broad-band colours (U — V), 
(U-B),(B-V), (V-R) and (V-K), the bolometric cor- 
rections, the specific contributions of evolutionary phases to 
the bolometric and to the broad band luminosities, and the 
stellar mass-to-light ratios (see App. A, Tables Al, A2 for a 
complete list). 
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Figure 7. The time evolution of the relative contributions 
L\P l /L\ "rj, of stars in the various evolutionary stages to the in- 
tegrated bolometric light of a stellar population. 



Figure 9. The same as Fig. 8 for the Main Sequence. 
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Figure 8. The time evolution of effective temperatures T e in the 
evolutionary PMS stages. 



4.1 The time evolution of integrated colours 

The time evolution of synthetic broad-band colours is shown 
in Fig. (jjj 

(U — V), (U—B) and (B—V) colours vary smoothly with 
age, because they trace the turnoff effective temperature. 
The evolutionary trend of IR colours is instead controlled 
by the AGB and RGB phase transitions. Their impact is 
particularly evident in (V — K), rising from ~ 1.3 to ~ 3.2 
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Figure 10. The time evolution of synthetic broad-band colours 
for a Salpeter IMF. Dashed lines represent models in which the 
TP-AGB contribution is not taken into account. 



in a short time interval. Its early rise is due to the AGB 
phase transition starting at t ~ 250 Myr. The RGB phase 
transition takes place at t ~ 500 Myr, contributing to keep 
values high of (V — K). The modest decrease in (V — K) 
at 1 Gyr Ss t ^ 4 Gyr reflects the reduction of the TP- 
AGB contribution. The TP-AGB phase plays a major role 
in determining the IR colours at intermediate ages. This 
can be appreciated looking at models in which the TP-AGB 
contribution has been omitted (also shown in Fig. |lFJ). The 
impact of the two phase transitions on optical colours is 
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Figure 11. The evolutionary phase contributions to the total monochromatic luminosities, L^/L^ T , as functions of the SSP age, for a 
Salpeter IMF. 



much weaker as AGB and RGB stars radiate mostly in the 
near infrared. 

The evolutionary behaviour of integrated colours can be 
easily understood by looking at the luminosity contributions 
of evolutionary stages, L»/Lj| T) shown in Fig. jn]. The UBV 
emission mainly comes from MS stars at any age, except for 
t 60 Myr, when the major V contributors are He-burning 
supergiants of intermediate spectral type. The K band re- 
veals the AGB phase transition at intermediate ages, when 
the TP-AGB contribution rises from ~ 20 per cent to ~ 80 
per cent. The RGB phase transition is visible in all bands, 



but particularly in K, where the RGB contribution domi- 
nates for t > 3 Gyr. The IMF has no influence on integrated 
optical colours (Fig. |l2"| ), as the variations induced in U, B 
and V luminosities are nearly identical, leaving the corre- 
sponding colours unchanged. A small dependence is visible 
only in the infrared, at ages greater than ~ 1 Gyr, when a 
dwarf-dominated (a; = 2.5) SSP has colours redder than a 
Salpeter one, because of the contribution of low-mass stars. 
The Main Sequence contribution is the most sensitive to the 
IMF slope (see Fig. ^)). 
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Figure 12. The influence of IMF on integrated broad-band 
colours: a Salpeter (solid line), a dwarf-dominated (dotted line) 
and a giant-dominated (dashed line) IMF. 



Figure 14. The relation between the scale factor A and the bolo- 
metric luminosity of a SSP, A/L^ oi (panel a). Also shown are the 
bolometric correction factors (panels b,c,d), for three different 
choices of the IMF slope. Solid line refers to a Salpeter, dotted 
line to a dwarf-dominated (x = 2.5) and dashed line to a giant- 
dominated (x = 0.5) IMF. 
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Figure 13. The dependence of MS luminosity contributions on 
the IMF slope, from a dwarf-dominated (x=2.5) to a giant- 
dominated (a:=0.5) one; £=1.35 is the classical Salpeter and Scalo 
refers to a three-slope IMF, defined in the caption of Fig. 2. 



.4 



,dM 1 



dt 



(11) 



where B(t) = b(t)/L^ A (t) is the evolutionary flux per unit 
luminosity of the SSP (the 'specific evolutionary flux'), a 
quantity nicely independent of the IMF slope (RB86). The 
total bolometric luminosity L^ ol (t) is related to the total 
monochromatic ones L^(t), through the bolometric correc- 
tion factors Bc(t) 



L h T °\t) = B*(t)L*(t). 



(12) 



B^ are functions of metallicity, IMF slope and age and are 
the direct product of synthetic SSP models. So, through eq. 
(11) and (12), we get A as a function of the luminosity in a 
certain band. 

The ratio between the scale factor A and the bolomet- 
ric luminosity, and the bolometric correction factors B* are 
shown in Fig. [l4| all as functions of age and IMF slope. The 
Lx o1 /-^t ratios are fairly independent of the IMF slope, as 
far as the IMF is not too steep for masses below ~ 0.5 Mq 
(i.e., 1 + x < 2.5). This is a consequence of the very narrow 
mass range producing most of the luminosity at nearly any 
SSP ages (see sec. 2). As an example, by adopting a Salpeter 
IMF, we obtain A ~ 1.12 ■ L^° l ~ 2.5 • L? for a 15 Gyr SSP. 



4.2 The IMF scale factor 

The size of a stellar population is direclty connected to the 
scale factor A of the IMF. Thus it is useful to express A in 
terms of the integrated luminosity, which is also proportional 
to the size of a SSP and directly observable. To this aim, we 
solve eq. (4) for A, obtaining 



4.3 The stellar mass-to-light ratio 

At any age, an SSP is composed of 'living' stars and 'dead- 
remnants', according to the initial mass of stars. The stellar 
mass is given by the convolution of the present mass M(t) 
with the IMF. The present mass M(t) coincides with the 
initial mass (M(t) = M in ) for M in < M TO (t) while M(t) 
= Mr (the remnant mass) for M[ n > M-ro(i). The kind of 
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Figure 15. The evolution of stellar mass-to-light ratios for a 
Salpeter IMF. 



Figure 16. The dependence of stellar mass-to-light ratios on the 
IMF. 



remnant also depends on the initial stellar mass: here we 
adopt the following recipe (cf. Renzini & Ciotti 1993). For 
Mro(i) < Mi n < 8.5, the remnant is a white dwarf of mass 
Mr = 0.077 M in + 0.48; for 8.5 < M in < 40, stars leave a 1.4 
Mq neutron star; finally, stars with Mi n > 40 are assumed 
to leave massive black- holes of mass Mr = 0.5 M; n . 

Fig. |l5| shows the evolution of stellar mass-to- light ratios 
for the bolometric and U, B, V, R and K luminosities. The 
early time evolution of A'h/L is very fast, then, after an 
age of ~ 5 Gyr, it becomes mild. For example, M t /L^ al 
increases of a factor of ~ 15 in the first 1 Gyr from the 
burst and only of a factor of ~ 3 from 5 Gyr to 15 Gyr. Fig. 
|li| shows the dependence of M»/L on the IMF, Mi n f being 
constant. At any age, the M*/L ratios of a dwarf-dominated 
population [x = 2.5) are much higher with respect to those 
of a population with a Salpeter or flatter IMF. For t J> 1 Gyr 
a population with a top-heavy IMF (x = 0.5) has M*/L 
ratios larger than those of a Salpeter SSP, because of the 
larger number of heavy remnants. 

The M t /L trend is easily understood by looking at 
the evolution of the total stellar mass and its components 
(Fig. |l7j). The stellar mass of a dwarf-dominated SSP is 
mainly composed of living stars at nearly any age because 
of the long evolutionary time-scales of low mass stars. A 
giant-dominated (x=0.5) SSP, on the contrary, has a larger 
amount of 'remnants' from massive stars. 

We refer to a forthcoming paper for an analysis of the 
dependence of M*/L on the SSP chemical composition. 
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Figure 17. The evolution of the total stellar mass, the mass in 
living stars and the white dwarf component (from the top to the 
bottom). Solid line is for a Salpeter IMF, dotted line for a dwarf 
dominated IMF (x = 2.5), dashcd-linc for a giant dominated IMF 
(x = 0.5) and long-dashed line for a Scalo multislope IMF. 




5 COMPARISONS WITH STAR CLUSTERS 

Synthetic broad-band colours and the luminosity contribu- 
tions of evolutionary phases are now compared with the 
corresponding quantities for Globular Clusters, the natural 
counterparts of SSPs models. Population size effects limit 



somewhat this comparison in the case of IR colours, as each 
cluster contains at most a few TP-AGB stars, hence stochas- 
tic fluctuations affect IR colours such as (V — K). 
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Figure 18. (U — B)-(B — V) synthetic diagram, compared with LMC and Galactic globular 
cluster colours. LMC GCs data (from van den Bergh 1981) are represented according to the 
SWB classification scheme (see the text), as indicated by the labels at the bottom left. NGC 
5927, NGC 6440, NGC 6624, Terzan 5, Pal 8 Galactic clusters (full triangles) and NGC 6528 
and NGC 6553 Bulge clusters (open triangles) data are from Harris (1996). Reddening values 
have been taken from the authors and the typical error bar is indicated. T96 and BC96 solar 
models are also plotted. 



5.1 Colour-colour diagrams 

Fig. |l8| shows our synthetic (U — B) versus (B — V) relation, 
compared with LMC GCs and Galactic clusters data. Obser- 
vational (B — V) colours are corrected for reddening follow- 
ing the authors prescriptions; for Galactic clusters (U — B) 
colours, the relation E(U - B) = 0.8E(B - V) has been 
adopted. Tantalo et al. (1996; hereafter T96) and Bruzual 
& Chariot (1996, private communication; hereafter BC96) 
models are also plotted. 

LMC clusters are grouped according to the SWB classi- 
fication scheme and represented with different symbols (in- 
dicated by the labels in Fig. |l| and |(]). Types 1 to 3 
correspond to young clusters (t 120 Myr, little fill trian- 
gles); types 3.5 to 5.5 to intermediate ones (250 Myr ^ t Si 
2 Gyr, fill dots) and types 6 to 7 to old ones (t <; 3 Gyr, 
open circles). 

The discrepancy between the models and LMC data 
for (B — V) <; 0.6 is due to the increasing difference in the 
metal content. In fact, as already noticed, SWB types > 6 
are much more metal poor (Z ~ 10 -3 — 10 -4 ) compared to 
earlier type clusters (Z ~ 0.01). This explains why they are 
bluer in spite of being older. 



In the range (B — V) Ss 0.6, a good agreement is evident 
and the various models behave in much the same way. Galac- 
tic cluster colours are consistent with old ages, in agreement 
with the standard picture. 

Fig. |l!^ shows (V — K ) vs. (U — B) compared with a sam- 
ple of LMC GCs with available infrared photometry (from 
Persson et al. 1983). 

As it is well known, intermediate age LMC clusters ex- 
hibit a wide range of (V — K) values, from ~ 1 to ~ 3.5. 
This is the result of the combined effect of AGB + RGB 
development, also visible in a (V — K) vs. (B — V) dia- 
gram (Fig. EOT) . The synthetic colours of our calibrated SSP 
model well agree with the observed ones, in this age range. 
T96 and BC96 models fail to reproduce the observed trend 
of IR colours at intermediate ages. Possible reasons will be 
discussed in the next paragraph. 



5.2 Comparing different EPS 

Analysing the discrepancies between different EPS models 
is not a simple job. Many factors contribute to produce 
different results, such as the stellar tracks used, the pho- 
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Figure 19. (V— K) — {U — B) two-colour diagram for the LMC GCs of Fig. ^ with available 
infrared photometry (from Persson et al. 1983); reddening values are taken from the authors 
and the typical error bar is indicated. As in Fig. [l8| clusters are represented according to the 
SWB parameter (see labels at the bottom-right). T96 models for two different metal content 
and Bruzual & Chariot (1996, private communication) solar models are also plotted. The 
lowest solid line represents our SSP model in which the TP-AGB phase has been omitted. 



tometric calibrations adopted to transform the theoretical 
Hertzsprung Russel diagram into observables and the syn- 
thesis method itself. In this section we limit ourselves to 
comment the discrepancy between synthetic (V—K) colours 
at intermediate ages, as shown in Fig. |l| and [U The com- 
parison is made with T96 and BC96. 

Fig. shows various synthetic broad-band colours as 
functions of age. The infrared indices of the three EPS mod- 
els are quite different. The impact of AGB phase-transition 
at intermediate ages is visible only in our (V — K) model. 
T96 anticipate its occurence, and their (V — K) shows a 
jump at t ~ 10 8 yr. This is in contradiction with the obser- 
vations. For example, the LMC cluster NGC 1866 (t ~ 10 8 
yr) does not have a developed AGB: according to FMB90, 
its AGB bolometric contribution is only the 6 per cent of 
the total. 

T96 adopt the analytical prescriptions for the AGB 
phase as in Bressan, Chiosi & Fagotto 1994 (hereafter 
BCF94). These are based on the core mass-luminosity rela- 
tion from Boothroyd & Sackmann (1988) and do not taken 
into account the results of Blocker & Schonberner 1991 (see 
3.2), as emphasized in BCF94 (see BCF94 for more details). 

BC96 adopt a semi-empirical receipt (described in 



Chariot & Bruzual 1991) to determine the position of TP- 
AGB stars on the theoretical CMD diagram. In their Fig. 
4, the authors compare the predicted fractional luminosity 
contributions of the AGB with FMB90 data, as we do in Fig. 
|H and a satisfactory agreement is evident. However, the in- 
tegrated (V — K) of BC96 does not reach the high values 
observed in the intermediate-age LMC GCs. This is proba- 
bly due to the energetic associated with C and M type stars 
(see Table 2). Furthermore, Chariot, Worthey & Bressan 
(1996) underline the lack of carbon star spectra in Bruzual 
& Chariot (1993) (and BC96) computations. The impact, 
on IR colours, of the amount of fuel burned by C or M type 
stars will be discussed in sec. 5.4.2. 

At the end, it is clear that the major source of discrep- 
ancy between all models is the treatment of the AGB phase, 
as it controls the SSP infrared emission at intermediate ages 
(cf. Fig. 0). 

On the opposite, the evolution of (U — B) and (B — 
V) is very similar for all sets of models. This arises from 
optical indices mainly depending on turn-off stars (cf. Fig. 
|l l|). The modest differences arise from the different input 
stellar tracks. 
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Figure 20. (V — K) — (B — V) two-colour diagram for the same LMC clusters and models 
as in Fig. 
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Figure 21. A comparison between different broad band colour 
models: this work (solid line); T96 (dotted line); BC96 (dashed 
line). All models assume a solar chemical composition, a Salpctcr 
IMF, and a lower mass cutoff at 0.1 Mq. 



5.3 The luminosity contributions of evolutionary 
phases 

A fundamental and global test for SSP models is the compar- 
ison of the synthetic luminosity contributions of the various 
evolutionary stages with the same quantities as measured in 
stellar clusters (cf. Renzini & Fusi Pecci 1988) . We show the 
results of two tests of this kind. 

The first test is performed on an intermediate age LMC 
GC sample (Ferraro et al. 1995; hereafter F95), selected with 
the aim of exploring the effect of AGB and RGB phase- 
transitions on integrated luminosities. The observed AGB 
and RGB contributions to the bolometric and K luminosity 
are compared with the synthetic quantities. 

As emphasized by the authors, their cluster sample is 
rather limited (12 clusters in total) and the available star 
sample in each cluster is subjected to statistical fluctua- 
tions, especially in the case of AGB stars. AGB stars, when- 
ever present, dominate the integrated cluster light, actually 
driving the IR magnitudes and colours, despite their small 
numbers. Anyway, the very short lifetimes involved (~ 3.5 
Myr for the whole AGB, Iben & Renzini 1983) make highly 
fluctuating the AGB luminosity contributions from one clus- 
ter to another of nearly the same age. To limit this effect, 
clusters have further been grouped in age bins, following 
the procedure described in sec. 3.2. The LMC cluster NGC 
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Figure 22. The AGB contribution to the total K luminosity 
compared with LMC GCs data from F95. The solid line represents 
the contribution of TP- AGB, while the dotted one refers to the 
whole AGB (i.e, E-AGB plus TP-AGB). Error bars indicate the 
r.m.s. (see the text). 



Figure 23. The RGB contribution to the total K luminosity 
compared to LMC GCs data from F95. 



1987 (s=35) was not taken into account, because the AGB 
contribution is not properly defined in F95. 

The cluster ages are derived from their s-parameter ac- 
cording to the relation logt = 0.079s + 6.05 (Elson k. Fall 
1988), as in F95. This relation is based on stellar evolution- 
ary models by VandenBerg (1985), which adopt canonical 
convection (no overshooting). Our SSP models are based on 
a different set of stellar models (CCS92), which also adopt 
canonical convection. The two age-turnoff luminosity rela- 
tions are very similar (to within ~ 10 per cent) and therefore 
the Elson & Fall s-age relation is consistent with our SSP 
models. For an s-age relation from overshooting models see 
BCF94. 

The results for the AGB phase are shown in Fig. |22[ 
The solide line represents the synthetic TP-AGB contribu- 
tion to the total if-light, and the dotted line refers to the 
whole AGB (that is, E-AGB plus TP-AGB). The error bars 
indicate the r.m.s. 

The comparison shown in Fig. ^ suggests an agree- 
ment between the observed AGB phase transition and the 
synthetic luminosity contribution of the AGB phase, that, 
we recall, has been calibrated on the LMC GCs, as described 
in sec. 3.2. 

Fig. H and |^ show the RGB contribution to the K 
and bolometric luminosities, compared with the same quan- 
tities from F95. Our model simulates the onset and early 
development of RGB phase transition very nicely. 

For old ages (t ~ 15 Gyr), the test is made with the 
~ Z e Bulge clusters NGC 6528 and NGC 6553, recently 
observed with HST (Ortolani et al. 1995). 

The empirical contributions of evolutionary stages to 
the V-luminosity (Lanteri Cravet et al. 1997 and Lanteri 
Cravet, private communication) are compared with the cor- 
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Figure 24. The RGB bolometric contribution compared to F95 
data. The synthetic £|^q b / ' L^° l takes into account that F95 sam- 
ple is limited to 2.3 mag. fainter than the TIP of the RGB. 



responding theoretical quantities in Table 5. To allow a more 
significative comparison, the synthetic MS ^-luminosity has 
been truncated at the completeness limits of the samples 
(~ 1.29 mag. below the TO for NGC 6553, for which 
ViFm 53 = 21.5, Ortolani, private communication; ~ 2.95 mag. 
below the TO, for NGC 6528, for which V^ 8 = 23.6, Mar- 
coni, private communication) . The theoretical contributions 
in Table 5 are fairly insensitive to the adopted IMF (see Ta- 
ble 6). In fact, down to the Vn m values quoted above, the stel- 
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Table 5. The contributions of evolutionary phases to the total V-luminosity (for a 
Salpeter IMF) compared with the observational data of the Bulge clusters NGC 6528 
and NGC 6553 (Lanteri Cravet et al 1997 and lanteri Cravet private communications) 



Phase j 


NGC 6528 


10 Gyr 


15 Gyr 


NGC 6553 


10 Gyr 


15 Gyr 


MS + SGB 


0.492±0.006 


0.457 


0.483 


0.421 ±0.04 


0.402 


0.432 


RGB 


0.271 ±0.01 


0.274 


0.273 


0.320 ±0.01 


0.302 


0.300 


HB 


0.154 ±0.02 


0.181 


0.163 


0.180 ±0.01 


0.199 


0.179 


AGB 


0.083 ±0.02 


0.087 


0.081 


0.08 ±0.01 


0.096 


0.089 



Table 6. The IMF influence on evolutionary phase luminosity contributions of Table 6 

AV^g 1 = 2.95 AKj,™ = 1.29 

x = 2.5 x = 0.5 x = 2.5 x = 0.5 

Phase j 10 Gyr 15 Gyr 10 Gyr 15 Gyr 10 Gyr 15 Gyr 10 Gyr 15 Gyr 

MS + SGB 0.476 0.502 0.444 0.471 0.410 0.441 0.397 0.425 

RGB 0.265 0.263 0.281 0.280 0.298 0.295 0.305 0.303 

HB 0.175 0.157 0.185 0.167 0.197 0.176 0.201 0.181 

AGB 0.084 0.078 0.089 0.083 0.095 0.088 0.093 0.090 



lar counts sample a small mass-range (for the 15 Gyr SSP we 
have: 0.82 < M < 0.938 for AV^o = 1.29, 0.7 < M < 0.938 
for AVro ~ 2.95). The agreement between predicted and 
observed contributions is quite satisfactory, and an old age 
(~ 15 Gyr) seems to be favoured, in accordance with Or- 
tolani et al. (1995) conclusions. 

5.4 Two experiments 

Thanks to the modular structure of the code it is possible 
to explore the impact of changing one or more ingredients 
and/or assumptions, on integrated photometric properties. 
We show here the results of two tests: the simulation of an 
'intermediate' morphology for the Horizontal Branch at late 
ages, and the effect of a different contribution of M and 
carbon stars at intermediate ages. 

5.4.1 Effect of a Blue HB 

The location of the Helium-Burning (HB) stars of Globular 
Clusters on a HR diagram displays different morphologies, 
called 'red-HB' (RHB), 'blue-HB' (BHB) or 'intermediate 
HB' (IHB), depending on the observed extension towards 
blue. As well known, the HB morphology shows a corre- 
lation with cluster metallicity (the 'first parameter'), but 
exceptions are so numerous that a 'second parameter' is re- 
quired. The nature of this second parameter has not been 
firmly established, but commonly considered candidates in- 
clude age and/or the cluster dynamical state (see Fusi Pecci 
et al. 1993). On average, metal-poor clusters have blue HB 
stars, while metal-rich ones show red HB. However, the pres- 
ence of BHB stars in old, metal rich SSPs is theoretically 
plausible (see Greggio & Renzini 1990) and now observa- 
tionally confirmed for some metal rich clusters in the galac- 
tic Bulge (Rich et al. 1997). In this section we explore the 
effect on integrated colours of adding BHB stars on the 15 
Gyr SSP, while keeping constant the HB fuel consumption. 
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Figure 25. The effect on integrated colours of adding a BHB 
to the 15 Gyr SSP. Dotted lines are SSP models in which an 
intermediate HB morphology is simulated. 

In this way, we simulate an IHB, as the red component is 
maintained. For this test we modify the {T e ,g)ij matrix at 
ii=15 Gyr, j=HB, extending the temperature distribution 
to logT c ~ 4; a flat stellar distribution is assumed. The 
surface gravity for each effective temperature is derived as- 
suming a constant luminosity. 

The result is shown in Fig. As expected, the (U — V) 
colour is the most influenced by the introduction of blue HB 
stars. 
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Figure 26. The influence of carbon stars populating the TP- 
AGB phase on integrated (V — K). Solid and dotted lines are 
the models of Fig. hcj (including or not the TP-AGB contribution 
respectively) ; dashed line is a model in which the amount of TP- 
AGB fuel for C stars has been reduced by a factor of two. 

5.4-2 Varying the fraction of carbon stars 

The initial metal content Z of a stellar population plays a 
key role in determining the production and the characteris- 
tics of carbon stars during the AGB phase (RV81). By re- 
ducing Z, the abundance of oxygen in the envelope is lower, 
thus a lower amount of carbon has to be dredge-up in order 
to achieve C/O > 1 and create carbon stars. Therefore, a 
metal poor stellar population is expected to display much 
more carbon stars than a metal rich one of the same age. 
Passing from Z ~ 0.008 (appropriate for LMC) to Zq, the 
fraction of carbon stars is reduced by a factor of ~ 2. To 
estimate this effect on integrated colours, we have reduced 
by a factor of 2 the amount of fuel burned by C stars (see 
Table 2), assigning it to M stars. 

Fig. ^6] shows the result: the AGB phase-transition is 
still evident in the integrated (V — K) at intermediate ages, 
but the colour maximum value is reduced by ~ 0.4 mag. 



6 CONCLUSIONS 

In this paper we have used the Fuel Consumption Theorem 
approach to compute evolutionary synthesis models (EPS) 
for Simple Stellar Populations (SSPs). The computational 
code has a modular structure, in which the main ingredi- 
ents are kept as separate entries. This allows to investigate 
their impact on the results and to better assess the uncer- 
tainties. The TP-AGB contribution has been calibrated on 
LMC globular Clusters. A grid of SSP models, having a so- 
lar chemical composition and covering a wide age-range (30 
Myr - 15 Gyr), has been constructed and submitted to sev- 
eral observational tests. 

Synthetic broad-band colours, compared with LMC 



GCs data, show a fair agreement. In particular, the ob- 
served rise in (V — K) colour at intermediate ages is well 
reproduced, as a consequence of the adopted calibration of 
the fuel consumption along the AGB phase. Comparing our 
results with others in the literature (Tantalo et al. 1996 and 
Bruzual & Chariot (1996; private communication), we show 
that only a correct calibration of the AGB energetics al- 
lows to match simultaneously the observed optical and IR 
properties of intermediate-age clusters. 

A complete test of the luminosity contributions by the 
various synthetic evolutionary phases has been performed 
on the old Bulge clusters NGC 6553 and NGC 6528. The 
result is very satisfactory and an old age (t ~ 15 Gyr) for 
these clusters is favoured. 

Having so tested the procedure, these models along with 
similar ones for other chemical compositions now being com- 
puted, will be used for a variety of astrophysical applica- 
tions. 
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APPENDIX A: TABLES 

The SSP models listed in the following tables refer to the 
TP-AGB fuel calibration described in Sec. 3.2. 
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Table Al. Integrated colours and stellar mass-to-light ratios as functions of age and IMF slope 
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Table A2. A/L^ and the bolomctric correction factors B£, as functions of age and IMF slope 

AGE(Gyr) A/L T B c u B C B B c v B C R Bf 
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